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ABSTRACT 

We investigate a stationary pair production cascade in the outer magnetosphere of 
a spinning neutron star. The charge depletion due to global flows of charged particles, 
causes a large electric fleld along the magnetic fleld lines. Migratory electrons and/or 
positrons are accelerated by this fleld to radiate gamma-rays via curvature and inverse- 
Compton processes. Some of such gamma-rays collide with the X-rays to materialize as 
pairs in the gap. The replenished charges partially screen the electric fleld, which is self- 
consistently solved together with the energy distribution of particles and gamma-rays at 
each point along the fleld lines. By solving the set of Maxwell and Boltzmann equations, 
we demonstrate that an external injection of charged particles at nearly Goldreich- Julian 
rate does not quench the gap but shifts its position and that the particle energy distri- 
bution cannot be described by a power-law. The injected particles are accelerated in the 
gap and escape from it with large Lorentz factors. We show that such escaping particles 
migrating outside of the gap contribute signiflcantly to the gamma-ray luminosity for 
young pulsars and that the soft gamma-ray spectrum between 100 MeV and 3 GeV ob- 
served for the Vela pulsar can be explained by this component. We also discuss that the 
luminosity of the gamma-rays emitted by the escaping particles is naturally proportional 
to the square root of the spin-down luminosity. 



Subject headings: gamma-rays: observations - gamma-rays: theory - magnetic flelds - 
methods: numerical - pulsars: individual (Geminga pulsar, PSR B1055-52, PSR B1706- 
44, Vela pulsar) 
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1. Introduction 

Recent years have seen a renewal of interest in 
the tlicory of particle acceleration in pulsar magne- 
tosphcrcs, after the launch of the Compton Gamma- 
ray Observatory (e.g., for the Vela pulsar. Kanbach 
et al. 1994, Fierro et al. 1998; for PSR B1706- 
44, Thompson et al. 1996; for Geminga, Mayer- 
Hassclwandcr et al. 1994, Fierro ct al. 1998; for PSR 
B1055 52, Thompson et al. 1999). The modulation 
of the 7-ray light curves testifies to the particle accel- 
eration cither at the polar cap (Harding, Tadcmaru, 
& Esposito 1978; Daughcrty & Harding 1982, 1996; 
Sturncr, Dermer, & Michel 1995), or at the vacuum 
gaps in the outer magnetosphere (Cheng, Ho, & Ru- 
derman 1986a,b, hereafter CHR; Chiang & Romani 
1992, 1994; Romani & Yadigaroglu 1995; Higgins & 
Henriksen 1997, 1998). Both models predict that elec- 
trons and positrons are accelerated in a charge de- 
pletion region, a potential gap, by the electric field 
along the magnetic field lines to radiate high-energy 
7-rays via the curvature process. However, there is 
an important difference between these two models: 
An polar-gap accelerator releases very little angular 
momenta, while an outer-gap one could radiate them 
efficiently. In addition, three-dimensional outer-gap 
models commonly explain double-peak light curves 
with strong bridges observed for 7-ray pulsars. On 
these grounds, the purpose of the present paper is 
to explore further into the analysis of the outer-gap 
accelerator. 

In the CHR picture, the gap is assumed to be 
geometrically thin in the transfield direction on the 
poloidal plane in the sense D± <C W, where D± repre- 
sents the typical transfield thickness of the gap, while 
W does the width along the magnetic field lines. In 
this limit, the acceleration electric field is partially 
screened by the zero-potential walls separated with 
a small distance D±; as a result, the gap, which is 
assumed to be vacuum, extends from the null sur- 
face to (the vicinity of) the light cylinder. Here, the 
null surface is defined as the place on which the local 
Goldreich-Julian charge density 



vanishes, where CI refers to the angular frequency of 
the neutron star, the magnetic field component 
projected along the rotational axis, and c the velocity 
of light; s and z refer to the coordinates parallel and 
perpendicular, respectively, to the poloidal magnetic 



field. The star surface corresponds to s = 0; s in- 
creases outwardly along the field lines. The last-open 
field line corresponds to 2 = 0; z increases towards 
the magnetic axis (in the same hemisphere). 

If i?2 > holds in the starward side of the null 
surface, a positive acceleration field arises in the gap. 
The light cylinder is defined as the surface where the 
azimuthal velocity of a plasma would coincide with 
c if it corotated with the magnetosphere. Its radius 
from the rotational axis becomes the so-called 'light 
cylinder radius', w^c = c/Q. Particles are not al- 
lowed to migrate inwards beyond this surface because 
of the causality in special relativity. 

It should be noted that the null surface is not a 
special place for the gap electrodynamics in the sense 
that the plasmas are not completely charge-separated 
in general and that the particles freely pass through 
this surface inwards and outwards. Therefore, the 
gap inner boundary is located near to the null surface, 
not because particle injection is impossible across this 
surface (as previously discussed), but because the gap 
is vacuum and transversely thin. 

Then what occurs in the CHR picture if the gap 
becomes no longer vacuum? To consider this problem 
rigorously, we have to examine the Poisson equation 
for the electrostatic potential. In fact, as will be ex- 
plicitly demonstrated in the next section, the original 
vacuum solution obtained in the pioneering work by 
CHR cannot be applied to a non-vacuum CHR pic- 
ture. We are, therefore, motivated by the need to 
solve self-consistently the Poisson equation together 
with the Boltzmann equations for particles and 7- 
rays. Although the ultimate goal is to solve three- 
dimensional issues, a good place to start is to examine 
one-dimensional problems. In this context, Hirotani 
and Shibata (1999a, b, c; hereafter Papers I, II, HI) 
first solved the Boltzmann equations together with 
the Maxwell equations one-dimcnsionally along the 
field lines, extending the idea originally developed for 
black-hole magnetospheres by Beskin et al. (1992). 
In Paper I, II, and HI, they assumed that the gap is 
geometrically thick in the transfield direction in the 
sense D± ::g> W. 

There is one important finding in this second pic- 
ture: The gap position shifts if there is a particle in- 
jection across either of the boundaries (Hirotani & 
Shibata 2001, 2002a,b; hereafter Papers VII, VIII, 
IX). For example, when the injection rate across the 
outer (or inner) boundary becomes comparable to 
the typical Goldreich-Julian value, the gap is located 



2 DIFFICULTIES IN PREVIOUS OUTER-GAP MODELS 



3 



close to the neutron star surface (or to the light cylin- 
der). In other words, an outer gap is not quenched 
even when the injection rate of a completely charge- 
separated plasma across the boundaries approaches 
the typical Goldreich-Julian value. Thus, an outer 
gap can coexist with a polar-cap accelerator; this 
forms a striking contrast to the first, CHR picture. It 
is also found in the second picture that an outer gap is 
quenched if the created particle density within the gap 
exceeds several percent of the Goldreich-Julian value. 
That is, the discharge of created pairs is essential to 
screen the acceleration field. 

The purpose of this paper is to examine the sec- 
ond picture more closely. In all the previous works 
in the second picture, the particle energy distribution 
has bcK;n assumed to be mono-energetic in the sense 
that the particles attain the equilibrium Lorentz fac- 
tor at each point, in a balance between the electro- 
static acceleration and the radiation-reaction forces. 
In the present paper, we discard this assumption and 
explicitly consider the energy dependence of particles 
by solving the Boltzmann equations for positrons and 
electrons. We will demonstrate that the particle en- 
ergy distribution cannot be represented either by a 
power law or by the mono-energetic approximation. 
We will further show that a soft power-law spectrum 
is generally formed in 100 MeV-GeV energies as a 
result of the superposition of the curvature spectra 
emitted by particles migrating at different positions. 

In the next section, we describe the difBculties of 
electrodynamics found in the first picture. We then 
present the basic equations in § 3. and apply the the- 
ory to four 7-ray pulsars and compare the predictions 
with observations in § 4. In the final section, we dis- 
cuss the possibilities of the unification of our picture 
with the CHR picture, as well as the unification of 
the outer-gap and polar-cap models. 

2. Difficulties in Previous Outer-gap Models 

To elucidate the electrodynamic difficulties in pre- 
vious outer-gap models, we have to examine the Pois- 
son equation for the electrostatic potential, 4*. Since 
the azimuthal dimension is supposed to be large com- 
pared with D± in conventional outer-gap models, the 
Poisson equation is reduced to the following two- 
dimensional form on the poloidal plane 



The true charge density, po, is given by pc = p+ 
where /o+ and p_ represent the positronic and elec- 
tronic charge densities, respectively. In the CHR pic- 
ture, the magnetic field is supposed to have a single- 
signed curvature. As a result, the particle number 
density {p+ — P-)/e grows exponentially in z direc- 
tion. Because of this exponential growth of the parti- 
cle number density, it has been considered that most 
of the 7-rays are emitted from the higher altitudes 
(i.e., large z regions) in the gap. 

To explain the observed 7-ray luminosities with a 
small D±, one should assume that the created current 
density becomes the typical Goldreich-Julian value 
in the higher altitudes. That is, the conserved cur- 
rent density per magnetic flux tube with strength B, 
should satisfy 



cp+ 
B 



+ 



-cp- 
B 



2^ 



(3) 



= ^-K[Pe{s,z)- PG3{S,Z)]. (2) 



in the order of magnitude. However, such a copious 
pair production will screen the local acceleration field, 
—d^/ds, as the Poisson equation (2) indicates. 

This screening effect is particularly important near 
to the inner boundary. Without loss of any generality, 
we can assume —d'^/ds > 0. In this case, because of 
the discharge, only electrons exist at the inner bound- 
ary. (We may notice here that external particle injec- 
tions are not considered in the CHR picture.) Thus, 
we obtain Pc/B = p-/B ^ — f7/(27rc) in the order 
of magnitude. In the vicinity of the inner boundary, 
we can Fourier-analyze equation (2) in z direction to 
find out that the —d^'^/dz'^ term contributes only to 
reduce d{—d'if/ds)/ds, the s gradient of the acceler- 
ation field. Thus, a positive — pcj must cancel the 
negative pc to make the right-hand side be positive. 
That is, at the inner boundary, 

B 2ttc B B ^ 2ttc ^ ^ 

must be satisfied, so that the acceleration field may 
not change sign in the gap. It follows that the polar 
cap, where i?^ ~ i? holds, is the only place for the 
inner boundary of the 'outer' gap to be located, if the 
created particle number density in the gap is compa- 
rable to the typical Goldreich-Julian value (eq. [3]). 
Such a non-vacuum gap must extend from the polar 
cap (not from the null surface where pcj vanishes) to 
the light cylinder. We can therefore conclude that the 
original vacuum solution obtained by CHR cannot be 
applied to a non-vacuum CHR picture when there is 
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a sufficient pair pro(hK;tion that is needed to explain 
the observed 7-ray luminosity. 

To construct a self-consistent model, we have to 
solve equation (2) together with the Boltzmann equa- 
tions for particles and 7-rays. We formulate them in 
the next section. 



A^-l-(.s, r) and (s, F) represent the distribution func- 
tions of positrons and electrons, respectively, at po- 
sition s and Lorentz factor T. We evaluate the di- 
mensionless Goldreich- Julian charge density Bz/B™ 
in equation (9) at each s, by using the Newtonian 
dipole field (see eqs.[9]-[12] in Paper VIII for details). 



3. Basic Equations 
3.1. Poisson equation 

Neglecting relativistic effects, and assuming that 
typical transfield thickness, D_\_, is greater than or 
comparable to the longitudinal width, W, we can re- 
duce the Poisson equation for non-corotational poten- 
tial \I' into the one-dimensional form (Hirotani 2000b, 
Paper VI; see also § 2 in Michel 1974) 



Pe(s) 



^Bz{s) 



2-KC 



■ (5) 



As described at the end of § 3 in Paper VII, it is 
convenient to introduce the Debye scale length c/wp, 
where 



27rce ' 



me 



(6) 



and B™ is the magnetic field strength at the inner 
boundary; e designates the magnitude of the charge 
on an electron, trie the rest mass of an electron. Thus, 
we can introduce the following dimensionless coordi- 
nate variable: 



1.87 X lO^a ^ 



VlOSGy/ 



1/2 



(7) 



where VL2 = ^/{Wv'aA s"^). 

By using such dimensionless quantities, we can 
rewrite the equation (5) into 



dE^^ 
d^ 

where 



A2 



m 



n+dV 



Ti-dT 



e5'(s) 



(8) 

gin 

(9) 
(10) 



meC^ c/wp ' 

and the particle distribution functions are defined by 



_ 2nceN±{s,r)_ 



Q B{0 



(11) 



3.2. Particle Boltzmann Equations 

On the poloidal plane, particles migrate along the 
magnetic field lines. Therefore, at time t, the dis- 
tribution function N of particles obeys the following 
Boltzmann equation, 

^ + P . V7V + Fext • 1^ = S{t, s, p), (12) 
at rrieT op 

where p refers to the particle momentum, i^ext the 
external forces acting on particles, and S the colli- 
sion terms. In the present paper, i^cxt consists of the 
Lorentz and the curvature radiation reaction forces. 
Since the magnetic field is much less than the crit- 
ical value (4.41 x 10^"^ G), quantum effects can be 
neglected in the outer magnetosphere. As a result, 
curvature radiation takes place continuously and can 
be regarded as an external force acting on a particle. 
If we instead put the collision term associated with 
the curvature process in the right-hand side, the en- 
ergy transfer in each collision would be too small to 
be resolved by the energy grids (see description below 
eq. [21] for details). We take the 7-ray production rate 
due to curvature process into account consistently in 
the 7-ray Boltzmann equations. 

3.2.1. Stationary Boltzmann Equations 
Imposing a stationary condition 



9 ^ 9 „ 



(13) 



and neglecting the pitch-angle dependence, we can 
reduce equation (12) to (Appendix A) 



9n+ 
dn- 

'df' 



+ 



E\\ — 



E\i 



meC^Wp 



dn- 



= 5+(€,r), (14) 



-5_(c,r), (15) 



where J7p designates the angular velocity of particles 
due to £J X Bp drift. Since the drift motion due to the 
gradient and the curvature of Bp is negligible for typ- 
ical outer-gap parameters, flp coincides ft provided 
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that i?0 = and < = Bpw/zui^c hold, 

where w refers to the distance from the rotation axis. 
Radiation-reaction force due to the curvature radia- 
tion, Pcv /c, is given by 



Pcv 2e2r4 



(16) 



where (s) refers to the curvature radius of the mag- 
netic field line. For a justification of adopting the 
pure-curvature formula, see § 5.1. 

3.2.2. Collision terms 

We assume in this paper that 7-rays are either out- 
wardly or inwardly propagating along the local mag- 
netic field lines. For simplicity, we neglect the devia- 
tion of the 7-rays from the field lines due to magnetic 
curvature. In what follows, we assume that the soft 
photons are emitted from the neutron star and hence 
unidirectional at the gap. Then the cosine of the col- 
lision angle /z has a unique value /i+ or for out- 
wardly or inwardly propagating 7-rays, respectively; 
/i+ and are determined by the magnetic inclina- 
tion a; at each position ^. Under those assumptions, 
the source term can be expressed as 

^pS+{^,T) = - / de^r?7c(e^,r,/i+)n+(^,r) 

+ [ dr,7?rc(r„r,/i+)n+(e,ri) 

^ri>r 



+ 



577p(e^,r,/x_,_) S7?p(e^,r,/x_) 

-9+ H ^ 9- 



dT 



dT 



(17) 



a;p5_(^,r) = - / de^r/7c(e7,r,/x_)n_(C,r) 



L 



drir,!c{ri,T,i^_)n-{^,Ti) 



+ 



ri>r 



5'?p(e7'r,M+) , 9^p(e7,r,/x_) 

-9+ H 9- 



dT 



dT 



(18) 



where the dimensionless 7-ray distribution function 
g± are defined as 

9±(e,e,)^^G±(s,e,); (19) 

G± refers to the number of 7-rays per unit volume per 
unit dimensionless energy e-y — hv^/ rrieC^. We will ex- 
plicitly define the diff'erential redistribution function 



for the pair production, d'q-p/dT, by equation (30), 
after briefly describing the soft photon field in equa- 
tions (23)-(29). 

If we multiply dT on both sides of equation (17) or 
(18), the first (or the second) term in the right-hand 
side represents the rate of particles disappearing from 
(or appearing into) the energy interval TOeC^F and 
meC^(F -|- dT) due to inverse-Compton (IC) scatter- 
ings; the third (or the fourth) terms does the rate of 
particles created via collisions between the outgoing 
(or the ingoing) 7-rays and ambient soft photons. The 
positrons (or electrons) are supposed to collide with 
the soft photons at the same angle as the outwardly 
(or inwardly) propagating 7-rays; therefore, the same 
collision angle cos~^ yu+ (or cos^^ ^_) is used for both 
IC scatterings and pair production in equations (17) 
and (18). The collisions tend to be head-on (or tail- 
on) for inwardly (or outwardly) propagating 7-rays, 
as the gap approaches the star. 

The IC redistribution function rfl^{e^,T , p) rep- 
resents the probability that a particle with Lorentz 
factor r upscatters photons into energies between e-y 



and 



de^ per unit time when the collision angle 



is cos~^ jjL. On the other hand, ?7jQ(Fi, F, fj) describes 
the probability that a particle changes Lorentz factor 
from Ti to F in a scattering. Thus, energy conserva- 
tion gives 



(20) 



We introduce the following dimensionless IC redis- 
tribution functions, 

1 f^' 

Vic±A'') = — Vici'^^, bi,IJ-±)de^, (21) 

where the particle Lorentz factor is discretized as F = 
bi; the 7-ray energy is divided into —lo + l + lm = 
15 -I- Im bins as big = b-u = 10"-^^/* x 61, 6_i3 = 
10-"/-* X 61, 6_i2 = 10-13/4 X 61, . . ., 6_2 = 10-3/-* X 
bi, b-i = 10-2/4 X 61, bo = 10-1/4 X 61, bi = 10^, 
62 = bi+Ab,. ..,h = 6i-h(i-l)A5,. . ., 6i„ = bi+{lm- 
1)A6. We normally use /„,, = 256 and A6 1.5 x lO'l 
In general, r?7c±.»(0 can be defined by the soft photon 
flux dFs/dcs and the Klein- Nishina cross section ctkn 
as follows: 

1-/3m± 



??IC±,i(0 



de. 



dF, 



de^^ [' dnl 
Jbt-i de-y J-1 



da 



KN 



'^dej*dn* 



(22) 
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where jS = -^/l — l/r^ is virtually unity, the solid 
angle of upscattered photon, the asterisk denotes the 
quantities in the electron (or positron) rest frame. In 
the rest frame of a particle, a scattering always takes 
place well above the resonance energy. Thus, the clas- 
sical formula of the Klein-Nishina cross section can be 
applied for the present problem. The soft photon flux 
per unit dimensionless photon energy eg [s~^cm~^] is 
written as dPs/dcs- 

We assume in this paper that soft photons are sup- 
plied by the thermal radiation from the neutron star. 
If a soft emission component cannot be represented by 
a black body spectrum, we use the observed spectrum 
of dPg/dcs corrected for the interstellar absorption. If 
it can be represented by a blackbody component, we 
adopt 

dF^^^BjTl(rj^)\ (23) 



de. 



\ r / 



where r and tns refers to the distance from the star 

center and the neutron star radius, respectively; the 
Planck function [ergs s~^cm~^ster~^ergs~^] is related 
with the blackbody temperature {6 = kT/nieC^) as 



2(mec2) 



2\3 



c^h^ exp(es/(5) — 1 



(24) 



Substituting equations (23) into (22), and executing 
the integration over e^, we obtain (Appendix B) 



dec 



exp(es/(5j) - 1 J_i 



dx*fic, 



(25) 



where 



fic 



jz -\- (a;*)^ - 1], for Zi_i < z < Zi 



otherwise; 



1 



l-Fes*(l-a;*) 
hi hi 



Zi = 



r(l -h /3a;*/x*)es* 6;(1 /3a;*M*)es 



(26) 
(27) 

(28) 



[I* and eg* are related with and Cg by the Lorentz 
transformation as 

M* = f7^, es*=r(l-/3/x)eg. (29) 



Since ?77c±,i(0 defined at F = 6;, F is replaced with 
hi in equation (28). Moreover, Aj represents the ob- 
served emission area of the ith blackbody component 
normalized by 47rr^g. For example, if a soft pho- 
ton field consists of the whole neutron star surface 
emission with temperature Tg and a heated polar cap 
emission with area 7rr^ and temperature Th, we ob- 
tain Ax = 1.0, (5i = kTs/m^c^ and A2 = 7rr^/(47rr^s), 
82 = kT^/rrieC?. 

To define a dimensionless pair-production redistri- 
bution function, we consider the case when a 7-ray 
photon with energy nicC^hi materializes as a pair with 
Lorentz factor hi. Dividing the pair production rate 
[s~^] per Lorentz factor by Wp, we obtain 

Vp±A^) = —^ihi,hi,iJ,±) 



i~fi± r dF^dop 



(30) 



where the pair-production threshold energy is defined 

by 

9 1 

(31) 



eth 



1 — /U e-y 

The collision angles cos~^ /i± are determined by the 

geometry, in the same manner as ?/7c±-«' ^^'^ differ- 
ential cross section for pair production is given by 



dar, 



f3? 



CM 



dr 8 



1 + /^cm(2 ~ Mcm) ^/^cm(-*- /^cm)^ 



1 /^cmMcm (1 /^cmMcm)^ 



32) 



where ctt refers to the Thomson cross section and the 
center-of-mass quantities are defined as 



±- 



2r 



0f 



CM 



1 



(1 - fj,)ese^ 



(33) 



where 11 = H+ (or /z_) for -q^+^i (or -qp-^i). 

Let us now discretize the differential equations. 
Denoting as the quantity / evaluated at ^ = ^'^ 
and T = bi (or e-y = hi), the source terms (eq. [17] 
and [18]) can be represented as 



i<l 



i>l 



B 



+ E[^p+.'W^+f»+^p-.'W^-fi]- (34) 
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The dimcnsionlc'ss 7-ray distribution functions are de- 
fined above the energy meC^bi^ = 5.11 x 10 eV as 



9±,i 



2TTce 



(35) 



where i = Iq + \, . . . , —1, 0, 1, 2, ... , Im- On the other 
hand, the particle distribution functions are defined 
above T > 61 = 10^ as 



n± 



J, _ 2'Kce ( c ^f. 



(36) 



where I = 1,2, ... ,1^- Note that n± are normalized 
by B at each point , while g± are by -B(^™). 

Since the IC scatterings do not alter the particle 
number, the sum of the first and the second terms in 
the right-hand side of equation (34) over /, vanishes 
identically. Therefore, 



^(5+f-5_f)=0 



(37) 



holds, because the created number of positrons is 
identical to that of electrons. 

3.2.3. Current Conservation 

Let us, for the moment, consider the current con- 
servation problem. Since we are considering a pure 
curvature process (i.e., neglecting the pitch-angle de- 
pendence of n±), and since F 1 holds, all the 
positrons (or electrons) migrate with velocity +c (or 
— c). As a result, particles at ^ = ^'^"^ migrate to 
^ = ^'^ after time (^'^ — £,^~^)/c. Figure 1 shows this 
situation, describing the particle propagation along 
the characteristics, whose gradient is given by (see 
eq. 14) 



-=^ll 



(38) 



rUeC^UJp 

in the phase space (^,r). Adding the produced pairs 
between ^'^"^ and which are depicted by open cir- 
cles in figure 1, wc obtain 



k-l 



fc 
+ 1 



(39) 



for positrons, where D'^ = — Equation (39) 

holds even if IC scatterings, which do not change the 
total particle number, contribute. In the same man- 
ner, we obtain 



(40) 



for electrons. Adding both sides of those two equa- 
tions, and utilizing equation (37), we obtain 

5](n+f, + n_f,) = l](n+fri + n.^'). (41) 



On these grounds, we can conclude that the current 
density 

jtt = E("+'+"-') (42) 

is virtually kept constant for k (i.e., for position ^ = 
^'=), provided T > 1. 

We solve the advection-type differential equations (14) 
and (15) on the two-dimensional phase space (^, F) 
by the Cubic Interpolated Propagation (CIP) method 
(e.g., Yabe & Aoki 1991, Yabe, Xiao, & Utsumi 2001). 

3.3. Gamma-ray Boltzmann Equations 

In this subsection, we consider the 7-ray Boltz- 
mann equations. For simplicity, we assume that the 
density of outwardly (or inwardly) propagating 7-rays 
decreases (or increases) at the same rate with the 
magnetic field strength; this assumption gives a good 
estimate when W <C wlg holds. Under this assump- 
tion, the 7-ray Boltzmann equations become (e.g.. Pa- 
per VII) 

J dV [r/7c(e7' r, M±) + VCYie^,^)] N±{s, F), 



+ 



(43) 



where ricv is the curvature radiation rate [s~^] into 

the energy interval mcC^(e-y + rfe-^) by a particle mi- 
grating with Lorentz factor F. The first term in the 
right-hand side expresses the 7-ray absorption rate by 
pair production, while the second term does the emis- 
sion rate via IC scatterings and curvature radiation. 

Integrating equations (43) over between bi-i 
and bi. we obtain 



dg± 



d{\nB) 



9±,i T^Vp±,i{bi) ■ 9±,i 



d^ d^ 

± E hc±,i(0 + VCY±,i{l)] n±(^, bi), (44) 



where r]cv+,i{l) refers to the curvature radiation rate 
VGv{bi,bi) by a positron having Lorentz factor F = 
bi. In the same manner, r]cv-,i{l) does that by an 
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electron. Evaluating the right-hand side at ^ = 

we obtain the advanced distribution functions by the 

Euler method as follows: 

54 = 5±r+^'(^)P- (45) 

3.4. Boundary Conditions 

As we have described, the set of Maxwell and Botz- 
mann equations consist of equations (8), (9), (14), 
(15), and (45). The Poisson equation and the 7-ray 
Boltzmann equations are ordinary differential equa- 
tions, which can be straightforwardly solved by a sim- 
ple discretization. On the other hand, the hyperbolic- 
type partial differential equations (14) and (15) are 
solved by the CIP method. 

In this section, wc consider the boundary condi- 
tions to solve the set of Maxwell and Boltzmann equa- 
tions. At the inner (starward) boundary (^ = ^™), we 
impose (Paper VII) 

^ii(r) = o, v(r) = o, (46) 

) =0 (i = /o + l, . . . , -1, 0, 1, 2, . . . , Im), (47) 
f -)■'" for / = 1 

where the dimensionless positronic injection current 
across the inner boundary ^ = ^™ is denoted as j™. 
As we described at the end of § 3.2, current density 
is conserved along the magnetic flux tube if particles 
are accelerated well above the initial Lorentz factor. 
Thus, we obtain at ^ = ^™ 

E"-.'(n=itot-r. (49) 
I 

At the outer boundary {£_ = ^°"*), we impose 

Ei\{r')=0, (50) 

5i(^-t) = (i = Zo + 1, • • • , -1, 0, 1, 2, . . . , /„), 

(51) 

The current density created in the gap per unit flux 
tube can be expressed as 

jgap=jtot-r-j°"*. (53) 

We adopt igap, i'", and as the free parameters. 



We have totally 2l„i + 2{lm — Zo + 1) + 4 bound- 
ary conditions (46)-(52) for 2lm + 2{lm -Iq + 1)+2 
unknown functions n±i, g±i, ^, and Thus two 
extra boundary conditions must be compensated by 
making the positions of the boundaries and be 
free. The two free boundaries appear because E^^ = 
is imposed at both the boundaries and because jgap is 
externally imposed. In other words, the gap bound- 
aries (^™ and ^°"*) shift, if j"^ and/or varies. 

In practice, we give n+(^'",r) and n_(^"\r) by a 
continuous function of F near the lowest F bin, 61, to 
avoid numerical oscillation due to the step function 
represented by equations (48) and (52). Neverthe- 
less, detailed functional forms of the injected particle 
spectrum little affect the results. 

3.5. Gap Position vs. Particle Injection 

Let us briefly examine how the gap position changes 
as a function of j™ and For a detailed argument, 
see § 2.4 in Paper VII. Defining the particle spatial 
number density by 

drn±(C,r), (54) 
we can rewrite the Poisson equation (9) as 

^ = ^[%ap(0 + /null(0], (55) 

where the created charge density in the gap is defined 
as 

engap(0 = e {n+(0 - r " [n-iO " (56) 
and 

/„„u(0 = r-i°"* + ^; (57) 

the two-dimensional screening effect (i.e., the —t/j/A']_ 
term) is neglected in equation (55). The created 
charge density increases outward due to discharge and 
satisfies ngap(f ") = -jgap and ngap(^°"*) = jgap. We 
have rigap < and /nuii > in the inner part of the 
gap, while we have rigap > and /nuU < in the 
outer part. It follows from equation (55) that the 
created particles partially screen the original electric 
field. If igap = i'" - + B^(^")/B'" holds, we ob- 
tain dE\^\ld^ = at the inner boimdary; that is, the 
gap is quenched due to the overfilling of space charges. 

Then, how about the injected particles? Let us ne- 
glect the created particles by imposing jgap <C 1, and 
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consider the following three cases: 

• J™ = 

In this case, maximizes (i.e., /null vanishes) at the 
point where ~ holds. That is, the gap is located 
around the null surface. 

• j'" ~ 1, ~ 

In this case maximizes at the point where = 
—B holds. That is, the gap is located near to the light 
cylinder. 

• .f " ^ 0, - 1 

In this case maximizes at the point where B^ = B 
holds. That is, the gap is located near to the polar 
cap surface. 

In short, created particles in the gap can quench 
the gap, while injected particles consisting of a sin- 
gle charge only shift it. The gap center is located at 
the 'generalized null surface' on which /nuii vanishes, 
provided that the created particle density is negligible 
(i.e.,ifngap <C 1). Strictly speaking, injected particles 
emit 7-rays that can materialize as pairs; therefore, 
injected particles contribute for screening of in a 
sense. Nevertheless, is screened by the created 
particles discharged by itself, not by the originally 
injected ones that do not return. 

4. Application to Individual Pulsars 

To solve the set of the Maxwell and Boltzmann 
equations, we must specify the X-ray field, dFg/dcs, 
which is necessary to compute the pair-production re- 
distribution function (eq. [30]). In this paper, we use 
the X-ray fluxes and spectra observed for individual 
rotation-powered pulsars to specify dF^/des- In § 4.1, 
we summarize the observed properties of the X-ray 
field. Then we apply the theory to the Vela pulsar 
in § 4.2, to PSR B1706-44 in § 4.3, to the Geminga 
pulsar in § 4.4, and to PSR B1055-52 in § 4.5. We 
assume that the solid angle of the emitted 7-rays is 
1 ster throughout this paper. 

4.1. Input Soft Photon field 

We consider the photons emitted from the neu- 
tron star surface as the seed photons for (7-7) pair- 
production and IC scatterings. That is, we do not 
consider power-law X-ray components, because they 
are probably magnetospheric and beamed away from 
the accelerator. We evaluate the IR photon field, 
which is needed to compute the IC scattering rate, 
from the Rayleigh-Jeans tail of the surface (black- 



body) component. In table 1, we present the observed 
properties of the four 7-ray pulsars exhibiting surface 
X-ray components, in order of spin-down luminosity, 

-^spin • 

Vela (J0835-4513) From Chandra observations in 
0.25-8.0 keV, the spectrum of this pulsar is turned 
out to consist of two distinct component: A soft 
component and a hard, power-law component. The 
hard component is probably a magnetospheric ori- 
gin and hence will be beamed away from the gap. 
Thus, we consider only the soft component as the X- 
ray field illuminating the outer gap. This component 
can be modeled either as a blackbody spectrum with 
kT = 1.49MK and A^ = 0.045A»(d/0.25)2, or as a 
magnetic hydrogen atmosphere spectrum with effec- 
tive temperature kT = 0.68MK (Pavlov et al. 2001). 
Based on high-resolution Ca II and Na I absorption- 
line spectra toward 68 OB stars in the direction of the 
Vela supernova remnant, Cha, Sembach, and Danks 
(1999) determined the distance to be 250 ± 30 pc. 
B1706-44 (J1710-4432) Gotthelf, Halpern, Dodson 
(2002) reported a broad, single-peaked pulsed profile 
with pulsed fraction of 23%, using the High Resolu- 
tion Camera on-board the Chandra X-ray observa- 
tory. They fitted the spectroscopic data to find (at 
least) two components, e.g., a blackbody of kT = 
143 eV with A = 0.129^* and a power-law compo- 
nent with photon index of —2.0. We consider that 
the former component illuminates the gap efficiently 
and neglect the second one. We adopt d = 2.5 kpc as 
a compromise between the smaller dispersion-measure 
distance of 1.8 kpc based on the free electron model 
by Taylor and Cordes (1993) and the larger H I kine- 
matic distance of 2.4-3.2 kpc derived by Koribalski et 
al. (1995). 

Geminga (J0633-I-1746) The X-ray spectrum con- 
sists of two components: the soft surface blackbody 
with kTs = 50 eV and A^ = 0.22A4d/0.16f and 
a hard power law with a = —1.6 (Halpern & Wang 
1997). A parallax distance of 160pc was estimated 
from HST observations (Caraveo et al. 1996). 
B1055-52 (J1059-5237) Combining ROSAT and ASCA 
data, Greiveldinger et al. (1996) reported that the 
X-ray spectrum consists of two components: a soft 
blackbody with kT^ = 68eVandAs = 7.3A,(d/1.53)2 
and a hard blackbody with /cTh = 320 eV and Ah = 
2.3 X 10~'*(d/1.53)^. The distance is estimated to be 
1.53 kpc from dispersion measure (Taylor & Cordes 
1993). 



4 APPLICATION TO INDIVIDUAL PULSARS 10 



Table 1: Input thermal X-ray field 



pulsar 


Ig -^spin 


distance 


n 












model 


refs. 




ergs s~^ 


kpc 


rad s~^ 


G 


eV 




eV 








Vela 


36.84 


0.25 


70.4 


12.53 


128 


0.045 






blackbody 


1,2 












59 


1.000 






hydrogen atm. 


1 


B 1706-44 


36.53 


2.50 


61.3 


12.49 


143 


0.129 






blackbody 


3 


Geminga 


34.51 


0.16 


26.5 


12.21 


48 


0.160 






blackbody 


4,5 


B1055-52 


34.48 


1.53 


31.9 


12.03 


68 


7.300 


320 


10-3.64 


blackbody 


6 



1: Pavlov et al. 2001; 2: Ogelman et al. 1993; 3: Gotthelf, Halpern, Dodson 2002; 4: Halpern & 
Wang 1997; 5: Becker & Triimper 1996; 6: Greiveldinger et al. 1996; 



4.2. The Vela pulsar 

We first apply the theory to the Vela pulsar, using 
the X-ray field modeled by the blackbody spectrum in 
§ 4.2.1-4.2.3, and by the hydrogen atmosphere spec- 
trum in § 4.2.4. 

4-2.1. Acceleration Field and Characteristics 

Let us first consider the spatial distribution of the 
acceleration field, 

which is solved from the Poisson equation (9). We 
may notice here that s is related with ^ by equa- 
tion (7). For the Vela pulsar, a small value jgap = 
4.6 X 10~^ gives the best-fit spectrum (see § 4.2.3 for 
details). 

To compare the effects of particle injection, we 
present the distribution for the three cases of 
ji" = (sohd), 0.25 (dashed), and 0.50 (dash-dotted) 
in figure 2. The magnetic inclination is chosen to be 
Q!i = 75°. We adopt = throughout this paper, 
unless its value is explicitly specified. The abscissa de- 
notes the distance along the last-open field line (i.e., 
the s coordinate used in eqs. [1] and [2]) normalized 
by OTLC- 

As the solid line shows, E^^ is located around 
the null surface when there is no particle injection 
across either of the boundaries. Moreover, E^^ varies 
quadratically, because the Goldreich-Julian charge 
density deviates from zero nearly linearly near to the 
null surface. 

As the dashed and dash-dotted lines indicate, the 

gap shifts outwards as j'" increases. When j™ = 0.5 
for instance, the gap is located on the half way be- 



tween the null surface and the light cylinder. This 
result is consistent with Papers VII, VIII, and IX. 

In figure 3, we present the characteristics of partial 
differential equation (14) for positrons by solid lines, 
together with E\\ (^) when j'" = 0.25 and ai = 75° 
(i.e., the dashed line in fig. 2). We also superpose the 
equilibrium Lorentz factor that would be obtained if 
we assumed the balance between the curvature radi- 
ation reaction and the electrostatic acceleration, as 
the dotted line. It follows that the particles are not 
saturated at the equilibrium Lorentz factor in most 
portions of the gap. 

In the outer part of the gap where is decreas- 
ing, characteristics begin to concentrate; as a result, 
the energy distribution of outwardly propagating par- 
ticles forms a 'shock' in the Lorentz factor direction. 
However, the particle Lorentz factors do not match 
the equilibrium value (dotted line). For example, near 
the outer boundary, the particles have larger Lorentz 
factors compared with the equilibrium value, because 
the curvature cooling scale is longer than the gap 
width. Thus, we must discard the mono-cncrgctic ap- 
proximation that all the particles migrate at the equi- 
librium Lorentz factor as adopted in Papers I through 
IX. We instead have to solve the energy dependence 
of the particle distribution functions explicitly. 

The particles emit 7-rays not only inside of the 
gap but also outside of it, being decelerated by the 
curvature radiation-reaction force. The length scale 
of the deceleration is given by 



u.,,o.-.(^)"(-ii^)V) 
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Since the typical Lorentz factor is a few times of 10^, 
Icurv is typically much less than uj-^c- Therefore, the 
escaping particles lose their energy well inside of the 
light cylinder. 

4-2.2. Particle Energy Distribution 

As we have seen in the foregoing subsection, the 
distribution function of the outwardly propagating 
particles forms a 'shock' in the outer part of the gap. 
In figure 4, we present the energy distribution of par- 
ticles at several representative points along the field 
line. At the inner boundary (s — 0.184ci7Lc)j parti- 
cles are injected with Lorentz factors typically less 
than 4 x 10^ as indicated by the solid line. Par- 
ticles migrate along the characteristics in the phase 
space and gradually form a 'shock' as the dashed line 
(at s — 0.205ci7Lc) indicates, and attains maximum 
Lorentz factor at s = 0.228ci7LC as the dash-dotted 
line indicates. Then they begin to be decelerated 
gradually and escape from the gap with large Lorentz 
factors ^ 2.5 x 10^ (dotted line) at the outer bound- 
ary, s = s°"* = 0.241tnLc- 

Even though the 'shock' is captured by only a few 
grid points for the dash-dotted line in figure 4, the 
CIP scheme accurately satisfies the conservation of 
the total current density. 



jtot = n+(^)-hn_(^) 

— J ^ J ^ Jgap ~ J 



(60) 



For this case, jtot is accurately conserved at 0.25 
within 0.2% errors even at the 'shock'. 

4.2.3. Formation of Power-law Gamma-ray Spec- 
trum 

So far, we have seen that the outwardly propa- 
gating particles are not saturated at the equilibrium 
value and that such particles escape from the gap with 
sufficient Lorentz factors suffering subsequent cooling 
via curvature process. It seems, therefore, reasonable 
to suppose that a significant fraction of the 7-ray lu- 
minosity is emitted from such escaping particles. 

We present the 7-ray spectrum emitted from out- 
wardly propagating particles (i.e., positrons) for the 
Vela pulsar with ai — 75°, jgap = 4.6 x 10"*^, 
j™ = 0.25 in figure 5. The dashed line represents 
the 7-ray flux emitted within the gap, while the solid 
one includes that emitted outside of the gap by the 
escaping particles. Therefore, the difference between 
the solid and the dashed lines indicates the 7-ray flux 




Fig. 1. — Schematic picture of characteristics in the 
phase space (C,r), where the dimensionless coordinate 
^ is related with s by equation (7). Positrons (or elec- 
trons) move along these characteristics when > 
(or E'li < 0). The filled circles indicate the particles 
moved from ^ = ^fc-i ^ certain time inter- 

val. Created particles (open circles) are assumed to 
migrate immediately outwards after the birth. 



J0835-4513 




0.2 0.4 
distance along field line / pi_LC 

Fig. 2. — Spatial distribution of -d*/rf.s for j'" = 
(solid), 0.25 (dashed), and 0.5 (dash-dotted), for the 
Vela pulsar when ai = 75° and jgap = 4.6 x 10~^ 
and = 0. The abscissa designates the distance 
along the last-open field line normalized by the light 
cylinder radius. 
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J0835-4513 

inclination= 75deg 




0.2 0.22 0.24 

distance along field line / pi_LC 

Fig. 3. — Spatial distribution of —d'^/ds{s) (dashed) 
for the Vela pulsar when ai = 75°, j™ = 0.25, 
= 0, and jgap = 4.6 x 10"^ The soft photon field 
is approximated by a single blackbody component. 
The characteristics are also shown by solid lines. The 
equiUbrium Lorentz factor, which would be obtained 
if we balanced the curvature radiation-reaction force 
and the electrostatic force, is indicated by the dotted 
line. 



J0835-4513 



distance=0.25 kpc 




2x10' 
Lorentz factor 



3x10' 



Fig. 4. — Particle energy distribution at several points 
along the magnetic field lines for the same case as 
in figure 3. Initial spectrum (solid line) evolves to 
dashed, dash-dotted, and dotted lines, as particles 
propagate outwards. 



emitted by the particles migrating outside of the gap. 
For comparison, we plot the phase-averaged EGRET 
spectrum, which is approximated by a power law with 
a photon index —1.7 (Kanbach et al. 1994). 

It follows from the figure that the 7-ray spectrum 
in 100 MeV-GeV energies can be explained by the 
curvature radiation omitted by the escaping parti- 
cles. We adjusted the transfield thickness as Dx = 
0.16k7lc = 2.8W^ so that the observed flux may be 
explained. The luminosity of the 7-rays emitted out- 
side of the gap contribute 48% of the total luminosity 
5.08 X lO^^ergs s'^ between 100 MeV and 20 GeV. 
In another word, we do not have to assume a power- 
law energy distribution for particles (as assumed in 
some of the previous outer-gap models) to explain 
the power-law 7-ray spectrum for the Vela pulsar. 
This conclusion is natural, because a power-law en- 
ergy distribution of particles will not be achieved by 
an electrostatic acceleration, and because magneto- 
hydrodynamic shocks (i.e., real shocks) will not be 
formed in the accelerator. 

Because the X-ray field is dense for this young pul- 
sar, the pair-production mean free path, and hence 
the gap width becomes small (for details, see Hirotani 
& Okamoto 1998, Hirotani 2000a, Paper IV; Hirotani 
2001, Paper V). As a result, the potential drop in 
the gap 2.24 x lO^^ y is only 0.81 % of the electro- 
motive force (EMF) exerted on the spinning neutron 
star surface ^ /i/tULc = 2.79 X 10^5 V. Nevertheless, 
this potential drop is enough to accelerate particles 
into high Lorentz factors, lO'^'^. 

4-2. 4- Hydrogen Atmosphere Model 

We next briefly examine the 7-ray spectrum when 
the soft photon field is modeled as the hydrogen at- 
mosphere spectrum (Pavlov et al. 2001), rather than 
the blackbody spectrum. In figure 6, we present the 
resultant spectrum for the same parameter set as in 
figure 5, except for the fitting model for the observed 
X-ray field. We fix s™ instead of jgap for comparison 
purpose. It follows that IC flux increases 3 times while 
igap does only 2 times. This is because the hydrogen 
atmosphere model gives much low-energy (IR) pho- 
tons compared with the blackbody model, as can be 
understood from its small effective temperature and 
large emitting area. Nevertheless, the obtained GeV 
spectrum little differ from the blackbody model. This 
is due to the 'negative feedback effect' of the gap elec- 
trodynamics (Paper V), which ensures the existence 
of solutions in a wide parameter range and suggests 
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Fig. 5. — Computed 7-ray spectrum for the Vela pul- 
sar for the same case as in figure 3. The dashed hne 
represents the emission from the gap, while the solid 
one includes the flux emitted by the particles being 
decelerated via curvature process outside of the gap. 
The soft photon field is modeled as a blackbody. 
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Fig. 6. — Same figure as fig 5, except that the soft 
photon field is modeled as a hydrogen atmosphere 
spectrum. 
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the dynamical stability of the solutions. 

4-2.5. Solutions in a Wide Parameter Space 

With both blackbody and the hydrogen atmo- 
sphere models, wc can explain the observed 7-ray 
spectrum in a wide parameter space 45° < ai < 75° 
and 0.125 < < 0.25, by appropriately choosing 
jgap and D±. With increasing jgap (<C 1), the vF^ 
[Jy Hz] peak energy increases because of the increased 
W, while the sub-GeV spectrum becomes hard be- 
cause of the significant 7-ray emission within the ex- 
tended gap (rather than outside of it) above GeV en- 
ergies. On the other hand, D± affects only the nor- 
malization of the 7-ray flux; the 7-ray luminosity is 
proportional to D"^. 

Let us first fix j™ at 0.25 and consider how the 
best-fit values of jgap and D± depend on aj. As we 
have seen, they are jgap = 4.6 x 10~^ and D± = 
0.155n7LC = 0.84s'" for ai = 75°, where the gap inner 
boundary is located at s = s™ = 0.1835tJ7Lc- How- 
ever, the ratio D±/s™ increases with decreasing ai 
and becomes 1.06 for 45° (cf. 0.84 for a; = 75°). 
From geometrical consideration, we conjecture that 
D± should not greatly exceed s™; we thus consider 
that a\ > 45° is appropriate for 7™ = 0.25. A large 
Q!i is preferable to obtain a small D±/zui^c. How- 
ever, since the radio pulsation shows a single peak, 
we consider that ai is not too close to 90°. On these 
grounds, we adopted ai = 75° as a moderate value in 
this section. 

Let us next fix ai at 75° and consider how the best- 
fit values of jgap and D± depend on j'". The ratio 
D^/tulc increases with decreasing 7'" and becomes 
0.94 for j'" = 0.125 (cf. 0.84 for j'" = 0.25). This is 
because the decreased flux of the outwardly migrating 
particles (due to decreased j™) must be compensated 
by a large D± to produce the same 7-ray flux. Thus, 
we consider j™ > 0.125 is appropriate for ai = 75°. 
A large is preferable to obtain a small D±/wi,c- 
However, for j™ > 0.25, the gap is so extended that a 
significant 7-rays are emitted above GeV within the 
gap; as a result, the sub-GeV spectrum becomes too 
hard. On these grounds, we adopted j™ = 0.25 as a 
moderate value. 

For 45° < ai < 75°, 0.125 < < 0.25, and ap- 
propriately chosen jgap and D±, TeV flux is always 
less than 3 x 10^° JyHz. Thus, one general point be- 
comes clear: TeV flux is unobservable with current 
ground-based telescopes, provided that the emission 
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solid angle is as largo as 1 stcr and that the surface 
thermal (not magnetospheric) X-rays are upscattered 
inside and outside of the gap. Since the magneto- 
spheric X-rays will be beamed away from the gap and 
their specific intensity is highly uncertain, we leave 
the problem of the upscatterings of magnetospheric 
(power-law) X-rays untouched. 

4.3. PSR B1706-44 

We next apply the theory to a Vela-type pulsar, 
PSR B1706-44. To consider as small as pos- 

sible, we adopt a large magnetic inclination 75°. We 
compare vFjj spectra for the three cases: j™ = 0.4, 
0.2, and 0.1. 

To examine how the sub-GeV spectrum depends 
on j'", we fix the uF^ peak at the observed value, 
^ 2 GcV, by adjusting jgap appropriately. For the 
solid, dashed, and dash-dotted lines in figure 7, we 
adopt jgap = 2.2 X 10"^, 1.8 x 10"^, and 1.5 x 10"^, 
respectively, and d = 2.5 kpc. Moreover, the perpen- 
dicular thickness is adjusted so that the predicted fiux 
may match the observed value (3.2 x 10^^ JyHz) at 
1.4 GeV; for the solid, dashed, and dash-dotted lines, 
they are Dj_ = 1.05roLc = 3.5s'", OMw^c = 4.1s'", 
0.49roLc = 5.3s™, respectively. The sub-GeV spec- 
trum becomes harder with increasing j'", because the 
ratio of the fiux emitted outside of the gap and that 
emitted within it decreases as the gap extends with 
increasing j™. However, as the solid line indicates, 
the obtained sub-GeV spectrum is still too soft to 
match the observation. For a smaller j™ (< 0.1), the 
sub-GeV spectrum becomes further soft. For a non- 
zero the inwardly shifted gap is shrunk to emit 
smaller 7-ray flux, which results in a further greater 
D±/s™. On these grounds, the predicted sub-GeV 
spectrum becomes too soft or the 7-ray flux becomes 
too small (i.e., D± becomes too large) for any param- 
eter set of j'" and if ai = 75° and d = 2.5 kpc. 

For a small inclination (ai < 75°), the gap is lo- 
cated relatively outside of the magnetosphere, be- 
cause the null surface crosses the last-open held line 
at large distances from the star. Because the pair- 
production mean-free path becomes large far from 
the star, the gap is extended for a small magnetic 
inclination. In such an extended gap, particles sat- 
urate at the equilibrium Lorentz factor at the outer 
part (Takata ot al. 2002) and emit most of the 7-rays 
around the central energy of curvature radiation. As 
a result, a hard sub-GeV spectrum can be expected; 
however, we have to assume a large D± that exceeds 



zui,c if = 2.5 kpc. 

Nevertheless, if d is much less than 2.5 kpc and 
a-i < 75°, we can explain the observed spectrum 
with moderate D±. For example, for d = 1 kpc, 
a; = 45°, j'" = 0.4, and jgap = 3.8 x lO'^, the gap 
exists in 0.51ti7LC < s < O.lAwi^c and the resultant 
spectrum (dotted line in figure 7) matches the ob- 
servation relatively well with a reasonable thickness, 
D_L = 0.80II7LC = l-5s'". 

On the other hand, for a large inclination (ai > 
75°), the sub-GeV spectrum becomes softer than a; = 
75° case for the same j"\ and jgap- Therefore, 
the spectrum will not match the observation whatever 
distance we may assume. 

We can alternatively consider a CHR-like outer 
gap. Assuming a small D± (say, 0.05), we find that 
the gap extends along the field lines due to the screen- 
ing effect of the zero-potential walls. In most portions 
of this extended gap, particles are nearly saturated at 
the equilibrium Lorentz factor. As a result, the sub- 
GeV spectrum becomes hard and match the obser- 
vation with appropriate peak energy around 2 GeV. 
However, in this case, the 7-ray flux becomes too 
small to match the observed value, unless we adopt 
an unrealistic distance (say, 300 pc). 

In short, the phase-averaged EGRET spectrum for 
this pulsar cannot be explained by our current model 
for any combinations of ai, j'", j°"*, jgap, and D±, if 
d = 2.5 kpc. Therefore, we suggest a small distance 
(e.g., 1 kpc) for this pulsar. 

4.4. The Geminga Pulsar 

Let us apply the theory to a cooling neutron star, 
the Geminga pulsar. For a small ai (e.g., 45°), the 
gap is so extended that the outer boundary exceeds 
the light cylinder. For a larger ai, on the other hand, 
not only the outer-gap emission, but also a polar-cap 
one could be in our line of sight. Since there has been 
no radio pulsation confirmed, we consider a moderate 
magnetic inclination ai ~ 60°. 

Since the soft photon field is less dense compared 
with young pulsars like Vela or B1706-44, the gap is 
extended along the field lines. For the set of param- 
eters i'" = 0.25, j°"* = 0, and jgap = 8.0 x 10-^ 
which gives the best-fit 7-ray spectrum, we obtain 
W ~ 0.59roLc- We present the spatial distribu- 
tion of obtained for this set of parameters as the 
dashed line in figure 8, as well as the characteris- 
tics (solid lines) and the equilibrium Lorentz factor 
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(dotted line). In the outer part of this extended 
gap, pGJ gradually increases with the distance, s. 
As a result, -E||(s) deviates from quadratic distribu- 
tion and decrease gradually in the outer part of the 
gap. Because of this extended feature, particles are 
nearly saturated at the equilibrium Lorentz factor. 
In another word, the mono-energetic approximation 
adopted in Papers I-IX is justified for this middle- 
aged pulsar. Particle distribution function forms a 
strong 'shock', which is captured only with one grid 
point, in 0.4 < s/k7lc < 0.45. As a result, jtot fluc- 
tuates a little as figure 9 indicates. Nevertheless, it 
returns to the 0.2525 level, which is 1 % greater than 
the value it should be (0.250), as the characteristics 
begin to be less concentrated beyond s = 0.45cc7lc- 

In figure 10, we present the resultant 7-ray spec- 
trum. Because of the nearly saturated motion of the 
particles, they lose most of their energy within the 
gap. As a result, 7-ray luminosity associated with 
the escaping particles (5.3 x lO^^ergs s~^), is negli- 
gibly small compare to that emitted within the gap 
(1.11 X lO^^ergs s~^), which is represented by the 
dashed line. 

It follows from the figure that the observed, soft 
spectrum between 200 MeV and 6 GeV can be ex- 
plained by the present one-dimensional model. It is 
interesting to compare this result with what obtained 
for the Vela pulsar (figs. 5 and 6). Between 100 MeV 
and 1 GeV energies, both spectra are formed by the 
superposition of the curvature radiation emitted by 
the particles having different energies at different po- 
sitions. The important difference is that the parti- 
cles are saturated at the equilibrium Lorentz factor 
in the gap for the Geminga pulsar, while they are 
nearly mono-energetic but only decelerated via cur- 
vature process outside of the gap for the Vela pul- 
sar. Because the particles are no longer accelerated 
outside of the gap, they emit 7-rays in lower ener- 
gies compared with those still being accelerated in the 
gap. As a result, the 7-ray spectrum for the Vela pul- 
sar becomes softer than that for the Geminga pulsar. 
Extending this consideration, we can predict that a 
7-ray spectrum below GeV is soft for a young pulsar 
and tends to become hard as the pulsar ages. 

4.5. PSR B1055-52 

Let us finally apply the present theory to another 
middle-aged pulsar, B1055-52. To obtain a large 7- 
ray flux for an appropriately chosen set of j™, jgap, 
and D±{< s™), we adopt a large magnetic inclination. 
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Fig. 7. — Computed 7-ray spectra for PSR B1706- 
44. The solid, dashed, dash-dotted lines represent 
the spectra for = 0.4, 0.2, 0.1, respectively, when 
d = 2.5 kpc and ai = 75°, while the dotted one for 
jin = 0.4 when d = 1.0 kpc and ai = 45°. 
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Fig. 8. — Spatial distribution of—d^/ ds (dashed line) 
for the Geminga pulsar when ai = 60°, j" = 0.25, 
jout _ g^jj^ j^^^ _ g Q ^ 10-5. Particles are satu- 
rated at the equilibrium Lorentz factor (dotted line) 
for this middle-aged pulsar. 
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Fig. 9. — Total current density, jtot) for the same case 
as in figure 8. Even though the particle distribution 

function forms a strong 'shock' in the Lorcntz factor 



direction in OAwhc < s < 0.45ti7LC (see fig. 
is conserved relatively accurately. 
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Fig. 10. Computed 7-ray spectrum for the Gcminga 
pulsar for the same case as in figure 8. The dashed 
and solid lines represent the same components as fig- 
ure 5. 



ai = 75°. 

Since the acceleration field and the particle energy 
distributions are similar to the Geminga pulsar, we 
present only the computed 7-ray spectra for this pul- 
sar in figure 11. The solid and dashed lines represent 
the spectra for the distances 0.50 and 1.53 kpc, re- 
spectively. For d = 0.50 kpc, jgap = 4.5 x 10~^ is 
chosen so that the peak energy of curvature radiation 
may match the observed peak energy. In this case, 
the gap exists in 0.1450tj7lc < s < 0.3882ci7lc- For 
d = 1.53 kpc, Jgap = 4.0 X 10~^ is chosen and the gap 
exists in 0.1475k7lc < s < 0.3818ti7Lc. 

It follows from the figure that the solid line matches 
the observed flux with a reasonable transflcld thick- 
ness, D± = 0.61s'" = O.'dSW. On the other hand, for 
the dashed line, we have to choose D± = 0.25zui,c = 
1.7s'", which is unrealistic. The observed fluxes can- 
not be explained with acceptable gap width (e.g., 
D± < s™) no matter what we may adjust j™, 
and Jgap if d = 1.53 kpc. 

On these grounds, we conjecture that the distance 
1.5 ± 0.4 kpc determined from the dispersion measure 
(Taylor & Cordcs 1993) is too large and that a more 
closer distance, such as 500 pc derived from ROSAT 
data analysis (Ogelman & Finley 1993) or 700 pc esti- 
mated from a study of the extended nonthermal radio 
source around the pulsar (Combi, Romero, Azcarate 
1997), is plausible. This conclusion is consistent with 
the results obtained in Paper IX under the mono- 
energetic approximation that all the particles are sat- 
urated at the equilibrium Lorentz factor, which is jus- 
tified for middle-aged pulsars. 

5. Discussion 

In summary, we have quantitatively examined the 
stationary pair-production cascade in an outer mag- 
netosphere, by solving the set of Maxwell and Boltz- 
mann equations one-dimensionally along the mag- 
netic field. We revealed that an accelerator (or a 
potential gap) is quenched by the created pairs in 
the gap but is not quenched by the injected particles 
from outside of the gap, and that the gap position 
shifts as a function of the injected particle fluxes: If 
the injection rate across the inner (or outer) bound- 
ary approaches the typical Goldreich- Julian value, the 
gap is located near to the light cylinder (or the star 
surface). It should be emphasized that the particle en- 
ergy distribution is not represented by a power law, as 
assumed in some of previous outer-gap models. The 
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particles escape from the gap with siifBcient Lorentz 
factors and emit significant photons in 100 MeV- 
3 GeV energies via curvature radiation outside of the 
gap. The 7-ray spectnim inchiding this component 
explains the phase-averaged EGRET spectra for the 
Vela pulsar, the Geminga pulsar and PSR B1055-52 
between 100 McV and 6 GeV. TcV fluxes are im- 
observable with current ground-based telescopes for 
these pulsars. 

We show that synchro-curvature process can be ap- 
proximated by the pure- curvature one in the next sub- 
section. We then point out an implication to the 7-ray 
luminosity versus the spin-down luminosity in § 5.2, 
and discuss future extensions of the present method 
in §§ 5.3-5.5. 

5.1. Synchrotron vs. Curvature Processes 

Let us first examine the pitch angles of particles 
injected across the inner boundary. Imposing that 
the synchrotron cooling length scale, 

= rrnec^ ^inry 

is greater than the distance, ezulc, for the particles 
to migrate before arriving the outer gap, we obtain 



sinx < 7.2 X 10" 



eo.iT 



1/2 

' ^ B^o\ (62) 



where eo.i = e/0.1, T2 = T/IO^, ^3 = ^/(lO^rad s"^), 
and Bio = -B/(10^°G). Particles having initial pitch 
angles greater than this value will lose transverse mo- 
mentum to reduce the pitch angles below this value, 
while migrating the distance sw^c- Therefore, if the 
particles are supplied from the polar cap with T ^ 10^ 
for instance, we can safely state that the the pitch an- 
gles are less than 10^"'' rad. 

If a particle having such a small pitch angle is ac- 
celerated to r = 10^ in a weak magnetic field region 
B = 10^ G, the ratio between the synchrotron and 
the curvature radiation rate becomes 



Pc sm X 



-Pcv yVmcC^/eB 
= 7.7 X lO-^i^^^S^Pf 2 



Pc _ smx 
0.5ra7LC 10-^ 



m 



where Be = B/IO^G, Tr = T/IO^. It follows that the 
synchro-curvature radiation can be approximated by 



a pure-curvature one for the partic;les injected across 
the inner boundary. It is noteworthy that the den- 
sity of created particles is much less than that of in- 
jected ones (i.e., jgap ^ j'"); therefore, the synchro- 
curvature effects for the created particles, which have 
much greater pitch angles compared with the injected 
ones, are negligibly small for the gap electrodynamics 
as well for the resultant 7-ray spectrum. 

5.2. Gamma-ray vs. Spin-down Luminosities 

It should be noted that the emission from the es- 
caping particles attain typically 40% of the total 7- 
ray luminosity for young pulsars. Thus, it is worth 
mentioning its relationship with the spin-down lumi- 
nosity, 

igpin = -/Qfi oc (64) 

where the braking index n is related to the spin-down 
rate as 

Q = -A;0". (65) 

If the spin down is due to the magnetic dipole radia- 
tion, we obtain n = 3. 

The outwardly propagating particles escape from 
the gap with spatial number density (eq. [11]) 

n /?°"* 

iVout = (r+Jgap)^^, (66) 

where = i?(^°"*). Therefore, the energy carried 
by the escaping particles per unit time is given by 

Lesc = D^^cNoutTcscrricC^ , (67) 

where resc(~ 10'''^) refers to the Lorentz factor of 
escaping particles. Note that r^sc is essentially de- 
termined by the equilibrium Lorentz factor (dotted 
line in fig. 3) near the gap center. Since the equilib- 
rium Lorentz factor depends on the one-fourth power 
of i^ii , the variation of Fosc on pulsar parameters is 
small. We can approximate 5°"* as 



B" 



P 



(68) 



where r°"* refers to the distance of the outer boundary 
of the gap from the star center. Let us assume that 
the position of the gap with respect to the light cylin- 
der radius, r°"'/tz7LC, does not change as the pulsar 
evolves; this situation can be realized if j™ — is 
unchanged. Evaluating _B at r = 0.5ti7LC, we obtain 

2 



_ 4rescWeC „2 
-^esc — M 

7re 
0.5 



± 



OC 



-'spin 



(69) 
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where n = 3 is assumed in the second line. To de- 
rive this conclusion, it is essential that the particles 
are not saturated at the equilibrium Lorentz factor. 
Thus, the same discussion can be applied irrespective 
of the gap position or the detailed physical processes 
involved. For example, an analogous conclusion was 
derived for a polar-cap model by Harding, Muslimov, 
and Zhang (2002). It is, therefore, concluded that 
the observed relationship oc ispin'^'^ merely re- 
flects the fact that the particles are unsaturated in 
the gap and does not discriminate the gap position. 

Let us compare this result with what would be ex- 
pected in the CHR picture. Since the gap is extended 
significantly along the field lines in the CHR picture, 
particles are saturated at the equilibrium Lorentz fac- 
tor to lose most of their energies within the gap, 
rather than after escaping from it. We can therefore 
estimate the 7-ray luminosity as 

igap = (D^D^W) ■ TVout • Pcv (70) 

where refers to the azimuthal thickness of the 
gap, Pcv [ergs s~^ (particle)"^] represents the curva- 
ture radiation rate. Noting that the particle motion 
saturates at the equilibrium Lorentz factor satisfying 
Pgv/c = e{—(M>/ds), recalling that the acceleration 
field is given by -rf*/d,s « ilBD^/Ap^c in the CHR 
picture, and evaluating B a,t r = w^c, we obtain 

CX Lspin, (71) 

where n = 3 is assumed again in the second line. Even 
though the escaping particles little contribute to the 
7-ray luminosity in the CHR picture, it is worth men- 
tioning the work done by Crusius-Watzel and Lesch 
(2002), who accurately pointed out the importance of 
the escaping particles in the CHR picture, when we 
interpret oc i^pjj^ relation. 

As we have seen, the particles being no longer ac- 
celerated contribute for the 7-ray luminosity that is 
proportional to igpfn- Reminding that the particles 
migrate with larger Lorentz factors than the equilib- 
rium value in the outer part of the gap (see fig. 3), 
we can expect roughly half of the 7-ray luminosity is 
proportional to L^^^^ (mainly between 100 MeV and 
1 GeV) , and the rest of the half to Lspin (mainly above 
1 GeV). As a pulsar ages, its declined surface emis- 
sion results in a large pair-production mean free path, 
and hence W. Because |/9gj| oc becomes small in 



the outer part of such an extended gap, (^) deviates 
from quadratic distribution to decline gradually in the 
outer part (fig. 8). As a result, particles tend to be 
saturated at the equilibrium value. On these grounds, 
we can predict that the 7-ray luminosity tends to be 
proportional to Lgpin with age, deviating from L^^^^ 
dependence for young pulsars. 

In the present paper, we have examined the set of 
Maxwell and Boltzmann equations one-dimensionally 
both in the configuration and the momentum spaces 
(i.e., only ^ and F dependences are considered.) In 
the next three sections, we discuss the extension of 
the present method into higher dimensions 

5.3. Returning Particles 

If we consider the pitch-angle dependence of par- 
ticle distribution functions, we can compute the radia- 
tion spectrum with synchro-curvature formula (Cheng 
and Zhang 1996). Moreover, we can also consider 
the returning motion of particles inside and outside 
of the gap. The returning motion becomes particu- 
larly important when both signs of charge are injected 
across the boundary. For example, not only positrons 
but also electrons could be injected across the inner 
boundary from the polar-cap accelerator. If E\\ > 
for instance, the injected electrons return in the gap. 
This returning motion significantly affects the Pois- 
son equation, if their injection rate is a good fraction 
of the Goldreich- Julian value. 

It remains an unsettled issue whether an outer-gap 
accelerator resides on the field lines on which a polar- 
cap accelerator exists. To begin with, let us consider 
the case when the plasma flowing between the po- 
lar cap and the outer-gap accelerator is completely 
charge separated. Such a situation can be realized, 
for instance, if only positively charged particles are 
ejected outwardly from the polar cap while there is 
virtually no electrons ejected inwardly from the outer 
gap. Neglecting the pair production, current conser- 
vation law gives the charge density, per unit mag- 
netic flux tube as 



where v refers to the particle velocity along the field 
line, and jtot the conserved current density per mag- 
netic flux tube. At each point along the field line, 
Pe should match pcj- If the field line intersects the 
null surface, pe must vanish there; this obviously vi- 
olates the causality in special relativity. Therefore, a 
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stationary ejection of a completely charge-separated 
plasma from the polar cap can be realized only along 
the field lines between the magnetic axis and those 
intersecting the null surface at the light cylinder. On 
these grounds, it was argued that an outer-gap ac- 
celerator, which is formed close to the last-open field 
line, may not resides on the same field lines on which 
a polar-cap accelerator resides. This has been, in fact, 
the basic idea that an outer gap will not be quenched, 
bccaiise the particles ejected from the polar cap will 
fiow along the different field lines. This idea was wel- 
comed in outer-gap models, because a gap has been 
considered to be quenched if the external particle in- 
jection rate becomes comparable to the Goldreich- 
Julian value, which was proved to be incorrect in this 
paper. 

In general, however, the plasmas are not com- 
pletely charge separated and consist of both signs of 
charge (e.g., positrons and electrons). Such a situ- 
ation can be realized, for instance, if both charges 
are ejected outwardly from a polar-cap accelerator, or 
if positively charged particles are ejected outwardly 
from the polar cap while electrons are ejected in- 
wardly from the outer gap, or if there is a pair pro- 
duction between the two accelerators. In these cases, 
the velocities of both charges will be adjusted so that 
both the current conservation and pc = pgj are sat- 
isfied at each point along the field lines. Therefore, 
it seems likely that a polar-cap accelerator and an 
outer-gap accelerator reside on the same field lines. 

To examine if there is a stationary plasma flow 

between the polar cap and the outer gap, we must 
extend the present analysis into two dimensional mo- 
mentum space in the sense that the pitch-angle depen- 
dence of the particle distribution functions is taken 
into account in addition to the Lorentz factor de- 
pendence. For example, if both charges are ejected 
from the polar-cap accelerator, electrons will return 
in the outer gap, screening the original acceleration 
field in the gap, and violating the original balance 
of Pc = pcj outside of the gap. Because the re- 
turning motion of particles can be treated correctly 
if we consider the pitch-angle evolution of the dis- 
tribution functions, and because the pair production 
is already taken into account, our present method is 
ideally suited to investigate the plasma flows and 
distribution self-consistently inside and outside of the 
gap. 



5.4. Unification of Outer-gap Models 

In addition to the extension into a higher dimen- 
sional momentum space, it is also important to extend 
the present method into a two- or three-dimensional 
configuration space. In particular, determination of 
the perpendicular thickness, D± , is important to con- 
strain gap activities. There have been, in fact, some 
attempts to constrain D± in the CHR picture. Since 
—d^/ds is proportional to BD\, particles energies, 
and hence the 7-ray energies increase with increasing 
_Dx (for a fixed B). Zhang and Cheng (1997) con- 
strained £)_L, by considering the condition that the 7- 
rays cause photon-photon pair production in the gap. 
Subsequently, Cheng, Ruderman, and Zhang (2000) 
extended this idea into three-dimensional magneto- 
sphere and discussed phase-resolved 7-ray spectra for 
the Crab pulsar. In addition, Romani (1996) dis- 
cussed the evolution of the 7-ray emission efficiency 
and computed the phase-resolved spectra for the Vela 
pulsar, by assuming that BD\ declines as . How- 
ever, in these works, screening effects due to pair pro- 
duction has not been considered; thus, the obtained 
D^, as well as the assumed gap position along the 
magnetic field, are still uncertain. 

On the other hand, in our approach (picture), D± 
is not solved but only adjusted so that the 7-ray flux 
may match the observations. Therefore, the question 
we must consider next is to solve such geometrical 
and electrodynamical discrepancies between these two 
pictures. We can investigate this issue by extending 
the present method into higher spatial dimensions. 

5.5. Unification of Outer-gap and Polar-cap 
Models 

Electrodynamically speaking, the essential differ- 
ence between outer-gap and polar-cap accelerators is 
the value of the optical depth for pair production. 
In an outer-gap accelerator, pair production takes 
place via 7-7 collisions and its mean-free path is much 
greater than the light cylinder radius. Therefore, a 
pair production cascade takes place gradually in the 
gap. In such a gap, E\\ is screened out by the 'general- 
ized Goldreich-Julian charge density', /null (eq. [57]), 
which increases outwards if • B > 0. Since /nuii 
is negative (or positive) in the inner (or outer) part 
of the gap, there is a surface on which the right- 
hand side of equation (55) vanishes, as long as the 
injected current density is less than the Goldreich- 
Julian value. The gap is located around this 'gener- 
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alizcd null surface', which is explicitly defined in § 2 
of Paper VII. 

On the contrary, in a polar-cap accelerator, pair 
production takes place mainly via magnetic pair pro- 
duction, of which mean free path is much less than 
the star radius for a typical magnetic field strength 
{B ~ 10^^ G, say). As a result, a pair production 
avalanche takes place in a limited region, which is 
called as the 'pair formation front', in the gap (Faw- 
ley, Arons, & Sharlemann 1977; Harding & Muslimov 
1998, 2001, 2002; Shibata, Miyazaki, Takahara 1998, 
2002; Harding, Muslimov, Zhang 2002). In the pair 
formation front, a small portion of the particles re- 
turn to screen out iSy. Such a returning motion can 
be self-consistently solved together with S|| by our 
present method, if we implement the magnetic pair 
production and the resonant IC scattering redistri- 
bution functions in the source terms of the particles' 
and 7-rays' Boltzmann equations. We can execute the 
same advection-phase computation in CIP scheme; 
thus, all we have to do is to add these source terms in 
the non-advection-phase computation, which is not 
very difficult. Since analogous boundary conditions 
(e.g., = for a space-charge limited flow) will be 
applied, wc expect the present method is also appli- 
cable to a polar-cap accelerator. This is an issue to 
be examined in our subsequent papers. 
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A. Reduction of Particle Boltzmann Equa- 
tions 

In this appendix, we derive the Boltzmann equa- 
tions (14) and (15) from equation (12). To begin 
with, we consider the momentum derivative terms 
Fext ■ dpN in equation (12). Let us describe the mo- 
mentum vector p as 



P 



(Al) 



where (e||, e^, ej,) forms the orthonormal basis: e\\ 
and e_L are the unit vectors parallel and perpendicular 
to the local magnetic field on the poloidal plane, and 
is the azimuthal unit vector. Introducing polar 
coordinates, we express the components as 



P\\ = pcosx, 

= psinxcosx_L, 
P;, = psinxsinx_L, 



(A2) 



to obtain 
d 

d 

dp4, 



d 
dpz 



cosx 



d sinx d 
dp 



sni X cos X ± -Tr- 
ap 

sinx_L d 
p sin X dx± ' 
d 

smxsmxx-^ - 
op 



P dx 

d cosx cos x_L 9 



P dx 

cosx sin x_L 9 
P dx 



+ 



cos x± d 
psinx dx±' 



(A3) 



The external force acting on a particle can be ex- 
pressed as 



Next, we consider the first and the second terms in 
equation (12). Decoupling the toroidal velocity asso- 
ciated with the drift motion as 



V = V — VuflpEcf 



(A6) 



we obtain 
ON 



1- — =— • VA'^ = 

dt nieT dt 

dN ^ 



V- VN 



dN ^ dN _ dN i^dN 



ds 



dz 



. (A7) 



It is worth noting that the z and derivatives in 
equation (A7) denote the advection of A'' in the trans- 
field directions. In the present paper, however, we are 
interested in only s dependence in the configuration 
space. Thus, we further reduce the right-hand side 
of equation (A7), by utilizing the fact that particles 
are frozen-in. Assuming a stationary magnetosphere 
as represented in equation (13), and integrating equa- 
tion (AT) over the momentum space, we obtain 



/ 



dN 
'dt 



■vVN 



■ VN + 



vo d<b 



cfp, 

(A8) 

where Vp = i^^e^^ + VzE-z- 

Let us introduce the averaged particle velocity such 
that 



(*p) 



VpNd^p 



I 



NcPp 



{V4,) 



i4,Nd^p 



I 



(A9) 



Nd^p 



Then noting (Bekenstein & Oron 1978) 

.l(*P 



-B 



P5 



/,-, \ - I K^p)l p 



(AlO) 



-Fext = q 



ds 



X Br 



Pcv V 



, (A4) 



where v = p/{Tme), and q designates the charge on 
the particle. Substituting equations (A2) into equa- 
tion (A4), and using equation (A3), we obtain 



ext 



dN _ \ f d-^\ _ Pcv] dN 

dp 9 cosx y j c dp 

-q (-^\ smxdN _ qBpdl^ 
\ ds ) p dx P dx±' 



The first term on the right-hand side describes the 
energy dependence of N. 



where plus (or minus) sign is chosen for outwardly (or 
inwardly) propagating particles, we obtain 



/ 



dN 
'dt 



V ■ VN 



d-^p = ±\vp\B-V 



/ Nd'p 



(All) 

where V • B = is used. We can neglect the az- 
imuthal derivative in equation (All), if the azimuthal 
thickness is large compared with those on the poloidal 

plane, or if is small compared with Bp. Under 
these assumptions, equation (All) reduces to 



dN 
'dt 



+ vVN 



d^p 



±K\Bp^^ 



d f N 



d^p 
(A12) 
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On these grounds, we neglect z and (p dependence 
of the distribution functions and approximate equa- 
tion (A7) as 



9N 



Neglecting the x± derivative in equation (A5), and 
assuming l^p] = ccosx, we finally obtain 



o 9 fN±\ 



, Pcv 
±cosx-e(--j- — 



ds J p dx 



S{s,p,x), 



dN± 
dp 

(A14) 



where (or N-) denotes positronic (or electronic) 
distribution function, and S refers to the source 
term averaged in a gyration. If we neglect % de- 
pendence, we obtain equation (14) and (15), where 
p = meC\/r2 — 1 « mecT is used. (It may be worth 
noting that x is different from the pitch angle defined 
in the synchro-curvature formula derived by Cheng 
and Zhang 1996, because the particles azimuthally 
drift, while the synchro-curvature formula was de- 
rived when the guiding center moves along the mag- 
netic field.) 



J 1059-5237 




photon energy [eV] 

Fig. 11. — Computed 7-ray spectrum for B1055-52. 
Escaping particles little contribute for the luminosity. 
The solid and dashed lines correspond to = 0.50 kpc 
and d = 1.53 kpc, respectively. Parameters are com- 
monly set as DjL = 0.09tULc, i'" = 0.25, = 0. 



B. Inverse-Compton Scattering Redistribu- 
tion Function 

In this appendix, we derive the redistribution func- 
tion for the IC scatterings. The number of photons 
upscattered by a single positron or electron into the 
energy interval E* and E* + dE* in time interval dt* 
in the solid angle d^*{OZ, (j)!) is given by 



dN^ = dt* / dE, 



.dFl 
dE: 



where the asterisk denotes that the quantity is mea- 
sured in the positron (or electron) rest frame. The 
incident photon flux per unit incident photon energy 
[photons s~^cm~^erg~^] is given by 



dFl 
dE: 



1 

E* 



i:{E:,n:)cose*dfit, 



(B2) 



-iQj-g ij jg ii^Q specific in- 



where I* [ergs s~ cm~ stcr^ 
tensity, Q* the angle between the photon momentum 
and the normal vector of a plane across which the 
incident flux is measured (fig. B12). 



scattered photon direction 
^ projected on (y, z) plane 



^ \ incident photon 
direction, which 
is on (x,y) plane 




17* =(cose|, 3ine|,0) 

Vy ^(cosBy, sinSYCOsCPy/ sinOysinCpY) 

*** ** * 

X =l7g*l7y= cosOg cosBy+sinSg sinOy sinCPy 



Fig. B13. — Definition of angles in the positron rest 
frame. 



B INVERSE-COMPTON SCATTERING REDISTRIBUTION FUNCTION 



24 




Fig. B12. — Scattering geometry in the observer's frame (non-primed) and the positron rest frame (primed). 



When the particle is moving with Lorentz factor 
r = 1/ yT^^^, we can relate the particle rest frame 
quantities (with asterisks) with the observer's frame 
(without asterisks) as 

E* = r(i - /3 cos es)E^, E* = r(i - /? cos e^)E^, 

(B3) 



r2(l-/3cos6is)2' r2 (1-/3 cos 61.^)2' 

(B4) 



cos 6*3-/3 . sm 

cost's = ^ TT^ smt/. 

^ 1 - /3 cos 



r(l -/3cos6is) 

(B5) 

cosOrfOs 

^°^Q'"- = rV-/3cosg3)^ - ^^'^ 

Moreover, we obtain the following Lorentz invariant, 



r* T 

-'s _ -'s 



(B7) 



In the observer's frame, the soft photon flux [ergs 

s~^cm~2erg~-'^] is given by 

= J h{E,,n,)cosedfl, = 7rB4T) (^)', 

' (B8) 

where Bs is the Planck function is defined by equa- 
tion (24). Since 

dCt* cos e*^=dn cos (B9) 

E* Ef, 



holds for r ^ tnS) we obtain 

dFl^dFl^ B,{T) /rNs\2 
dE* dE, ^ E^ \ r ) ' 



(BIO) 



as long as the blackbody emission comes from the 
whole surface of the neutron star. If the observed 

emission area of the ith blackbody component is 
AjAnr-^g, equation (BIO) is modified as 

'''' Kt^)'^E^.-^^(^.) (BID 



dE; 



The IC redistribution function is defined by 



Vic 



(e^,r,jUs) = j dfl^ 



dN^ 



dtd{E^ /meC^)dCl-y 



(B12) 



where /is = cos^?s- Noting that 



dtdE^dn^ = r^{i - 13 cos 0^)dt*dE*dn*, (bi3) 

and assuming that the specific intensity of the inci- 
dent photons are unidirectional (i.e., jig = jj.^ =constant), 
we obtain 



7?7c(e7>r,Ms) = 



1 - /^/i^ 



dn* 



1 + /3 cos ( 



/ 



^ dca de*d^* ' 



(B14) 



where 



dFs ^ 27i^meC^ /rNs\2 ^ 



de. 



h^c 



exp{es/Sj) - 1' 
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. (B15) 



In the particle rest frame, the Klein-Nishina cross sec- 
tion is given by 



da 



KN 



de*(m* 



16^ ^ 



*\2 



+ {X*) 



X 5 



l + e*(l-x*) 



where 



TP* 

to , 



(B16) 



(B17) 



rUeC^ ' meC^ 

the quantity x* refers to the angle between the inci- 
dent and the scattered photon momenta (fig. B13), 
and can be expressed as 



X = cos 0g cos 0^ + sm t's sm 0^ cos cp^ . 



(B18) 



We define 5** to be the plane containing both the 
incident and scattered photon momenta. Then a* 
means the angle between the normal of S* and z axis. 
Denoting the normal vector as 



h = {K,hy, ^1-hl-hl) 



(B19) 



and imposing 



h v* = hx cos 01 + hy sin 0* = 0, 
h - V* = hx cos 0* + hy sin 0* cos (p* 

+ ^1-hl- hi sin 0*^ sin = 0(B20) 



we obtain 

cos a* = y^l — /i| — /i^ = 



/e 



\J /e^ + sin^ < 



:, (B21) 



where /e = sin 0*p cot 0* — cos 0* cos </)* . Solving equa- 
tion (B18) for cos<;i* and eliminating (^* in /e, we 
obtain 

6';-a;*cos6'* , ^ 

cos<^; = -^-^==4?^, (B22) 



sin 61*^1 - (a;*)2 



which leads to 



Since dfJ* = d(- cos = dx* da* holds, we 

can change the integration variables to obtain 



»?i'c(e7.r,Ms) = 



d.v" d( 



1 -f /3(sin^o*-\/l — (a;*)2cosa* -|- a;* cos^°*) 



X / de 



des dej*dn* 



(B24) 



Here, e* and e* in the right-hand side of equa- 
tion (B16) should be replaced with eg and by equa- 
tions (B3). 

We define the dimcnsionlcss IC redistribution func- 
tion by equation (21). Then, integrating over be- 
tween bi-i and bi, we obtain 



^7c,i(r) 



3ctx 



(l-/3/x°)/ de, 



dF, 



dx* 



/.27r 

/ da*fic{x*,a*,e,-n,,r), (B25) 
Jo 



where fic is defined by equation (26). Substituting 
dFs/dcs, and summing up the blackbody components, 
we obtain equation (25), where /j.^* and /U° are de- 
noted as /i* and fi in equation (25). 



cos 0* = sin 0* \/l — (x*)^ cos a* + x* cos 

w a;* cos 61* ifr»l. (B23) 



